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"by Milton R. Lauver, Jerry L. Hall, and Frank E. Belles 
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SUMMARY 

Gas temperature histories in a shock tube vere determined By an infrared 
monochromatic radiation pyrometer technique. A 10 percent carbon dioxide- 
argon mixture vas studied vith Both incident and reflected shocks. Measured 
temperatures of 1100° to 3300° K immediately "behind the incident shocks vere 
in good agreement vith the temperatures calculated "by one-dimensional vave 
theory from the shock velocity. Measured temperatures of 1900° to 3550° K "be- 
hind reflected shocks vere in good agreement for short reflection distances 
vith those calculated "by one -dimensional vave theory from the incident shock 
velocity. The relation "betveen the velocity attenuations and the temperature 
and pressure rises Behind incident shocks vas in general agreement vith strong 
shock theory. The change in gas temperature vith time Behind a reflected shock 
vas adequately calculated from the dissociation rate of carBon dioxide and the 
measured pressure history of the gas mixture. 


INTRODUCTION 

The shock tuBe has one of its prime applications in the study of high- 
temperature rate processes in gases. Commonly, these processes have a strong 
temperature dependence, and it is therefore necessary to knov the temperature 
of the experiment quite accurately. This need can Be met nicely so long as in- 
cident shock vaves are used to heat the gas, Because one-dimensional theory, 
plus the thermodynamic properties of the gas, can Be used to calculate the 
shocked gas temperature from the measured velocity of the vave. Very often, 
hovever, it is convenient or even necessary to study processes Behind reflected 
shocks, in order to take advantage of the higher temperatures and pressures that 
can Be produced vithout unduly straining the capacities of conventional shock 
tuBes. In that event, it is no longer so clear that the temperature can Be cal- 
culated accurately; interactions of the reflected shock vith the Boundary layer 
(ref. 1) and vith the pressure gradient created By attenuation of the incident 
shock (ref. 2) introduce complications. 

This problem has Been videly recognized, and measurements of reflected gas 
temperature have recently appeared in the literature. For the most part they 


were made by the line-reversal method (refs. 3 and 4), although in one case the 
rate of a chemical reaction was used as an indication of temperature (ref. 5). 
The purpose of the present work was to study reflected shock temperatures by a 
more direct method based on simultaneous measurements of infrared spectral radi- 
ance and emissivity. Among the advantages of the infrared method are the fol- 
lowing: (l) there is no need to introduce sodium salts or other sources of 

emission; (2) good time resolution is possible; and (3) a very wide temperature 
range can be covered, a range not limited by the available brightness tempera- 
ture of a light source as it is in the reversal method. 

This report gives temperatures measured behind reflected shocks in a gas 
that is typical of mixtures likely to be used in chemical studies, 10 percent 
carbon dioxide- 90 percent argon. The temperature range was 1900° to 3550° K. 

A series of temperatures measured at 50 to 200 microseconds after the reflected 
wave passed the observation station were extrapolated to zero time and compared 
with values calculated from the velocity of the incident wave. Two detailed 
temperature-time records are also presented; these show the effects of shock 
attenuation on both the incident and the reflected shock temperature. 


APPARATUS ARID EXPERIMENTAL PROCEDURE 
Basis of Emission-Absorption Pyrometry 

The infrared monochromatic radiation (IMRA) method of gas temperature 
measurement has been described elsewhere (refs. 6 and 7). It is based on the 
fact that the radiance of a hot gas depends on the gas temperature and on the 
number of gaseous entities (molecules, free radicals, or atoms) emitting radi- 
ation. Consequently, by measuring the absolute intensity of radiation (spec- 
tral radiance) and the relative number density of emitters (spectral emissiv- 
ity), the temperature of a gas can be calculated. Eigure 1 shows schematically 
the information recorded in a shock-tube experiment and the way in which it is 
reduced to a temperature measurement. 

The measurement of spectral emission is straightforward. The infrared de- 
tection system response to the hot gas emission V g (fig. 1) is compared with 
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Figure 1. - Basis of infrared monochromatic radiation pyrometer. 
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its response to a calibrated standard infrared source V g and converted to 

absolute units by means of calibration factors. The result is a spectral radi- 

-1 2 1 

ance g in terms of watts steradian centimeter” micron”* at a specific 

wavelength. The effect of the number of emitters is accounted for by absorp- 
tion spectroscopy at the same wavelength. The reduction in radiance of a 
chopped transmitted light beam from Vq to V volts determines the spectral 
absorptivity (1 - V /Vq ) , which is equal to the spectral emissivity e. 

Kirchoff’s radiation law relates the spectral radiance of a hot gas. to 
that of an ideal blackbody radiator. Planck’s law relates the spectral radi- 
ance of an ideal blackbody radiator to its temperature and emissivity. Combin- 
ing these laws results in the following relation between the spectral radiance 
of a gas N-v , its temperature T , and its spectral emissivity e: 

'v o § 



Thus, inasmuch as and C 2 are the known Planck radiation constants and 

g and e are measured at wavelength A, the gas temperature T^ may be 
calculated. 


Apparatus for Emission-Absorption I^rometry 

A schematic diagram of the shock tube and associated equipment is given 
in figure 2. A constant intensity of infrared radiation from a glower source 

G Glower 
B Internal standard 
C i 80-kcps chopper 

C 2 1. 1-kcps chopper 
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Figure 2. - Schematic diagram of infrared monochromatic radiation pyrometer. 
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is mechanically chopped at 80 kilocycles per second and sent through the sap- 
phire windows of the shock tube into a prism monochromator and a liquid- 
nitrogen-cooled indium antimonide detector. The shocked gases emit radiation, 
which passes out of the shock tube into the same radiation detection system. 

Its output is amplified "by a directly coupled circuit providing a measure of 
the radiance of the gas (the glower radiance is negligible in comparison to 
that of the gas). The glower signal is selected by a capacitor- coupled, 
80-kilocycle tuned amplifier of high gain. With these amplifiers, the detec- 
tion system uses frequency separation to distinguish between the emitted and 
transmitted light, although both are presented simultaneously to the one detec- 
tor, and sends the resulting signals to different oscilloscope beams for re- 
cording. 

The principle of temperature measurement used in the IMRA technique corre- 
sponds to that of the "Two-Path Method with Comparison Source" (ref. 8). An 
error analysis by Warshawsky (ref. 9) shows large errors for the latter as the 
gas temperature rises above the glower temperature. The crucial feature of the 
IMRA apparatus, which minimizes these errors, consists of chopping the glower 
radiation before it passes through the test gas. This permits the emitted and 
transmitted components of the radiation to be separated so that each may be am- 
plified independently by the amount required to obtain easily measurable oscil- 
lograms. Without chopping, the light from a source that is emitting more weak- 
ly than the gas would be swamped by the gas emission. 


Calibration of Pyrometer 

The IMRA apparatus was calibrated by determining the transmittance of the 
shock-tube window and the absolute radiance of the internal secondary- standard 
source (B in fig. 2), which was a tungsten-ribbon lamp. 

The window was found to pass 82 percent of 4.5-microns radiation falling 
on it. This was the wavelength used for the temperature measurements. Al- 
though the center of the asymmetric stretching band of carbon dioxide is at 
4.3 microns, the longer wavelength is much more desirable. At 4.5 microns, 
only the hot gas in the shock tube absorbs light, while the cool carbon dioxide 
in the room and in the shock-tube boundary layer is almost completely trans- 
parent. 

The internal standard was compared with a source that had been calibrated 
at the Bureau of Standards. It was found to emit 2.47 watts steradian" 1 

p i 

centimeter" micron” at 2.2 microns. The relative sensitivity of the optical 
system and detector at 4.5 and 2.2 microns was then determined hy measuring the 
radiation from a blackbody source of known temperature at the two wavelengths. 


With the foregoing calibration factors at hand, the absolute spectral 

radiance of the hot gas at 4. 5 microns is readily determined from the 

A, g 

following formula: 
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where 


Vg amplitude of emission signal (at 4.5 p), V 

V s amplitude of signal from internal standard source (at 2.2 p), V 

R sensitivity of the system at 4,5 relative to that at 2.2 p 

t w transmittance of shock-tube window at 4. 5 p 

g absolute spectral radiance of internal standard at 2.2 p 

This formula is valid provided the amplification and the slit width are 
the same when hoth Vg and V s are recorded. If they are not the same, the 

appropriate correction terms must be applied. 


Shock Tube 

The shock tube was of rectangular cross section, 37 by 74 millimeters. 

Two circular sapphire windows, 28 millimeters in diameter, were flush mounted 
opposite one another across the longer dimension. The midpoint of the windows 
was 3. 94 meters from the polyester plastic diaphragm that separated this driven 
section of the shock tube from the driver section. The midpoint of the windows 
was 179 millimeters from the downstream end wall of the driven section. This 
distance to the end wall could be reduced to 27 or 78 millimeters by means of 
close-fitting plugs. The diaphragms were pres sure -burst with helium. 


Timing, Pressure, and Recording Instrumentation 

Four thin- film resistance gages were mounted upstream of the windows, 
while one was mounted at the same axial position as the centerline of the win-’ 
dows. These gages marked the time of incident (and in some cases, reflected) 
shock arrival at the five positions, relative to the first, or trigger, posi- 
tion. The signals were displayed on one beam of a dual beam oscilloscope. The 
other beam displayed the output of a quartz pressure transducer. This trans- 
ducer was also located at a position corresponding to the centerline of the win- 
dows. A second dual-beam oscilloscope was triggered to display absorption and 
emission levels from the gas as measured by the IMRA apparatus. Timing pulses 
from a crystal- controlled secondary frequency standard were recorded on all 
four oscilloscope beams for each run. A calibration signal for the emission 
level interrupted about 1100 times per second for identification and a cali- 
bration signal for the pressure were also recorded each time. 

Figure 3 shows the oscillograph records of infrared absorption and emis- 
sion for a single experiment. The important features of this typical record 
are: (l) the large changes of transmitted and emitted light upon passage of 

the shock waves, (2) the relatively noise-free signals; and (3) the quick re- 
covery of the tuned amplifier in the absorption channel from the ringing in- 
duced by step changes, which permits meaningful readings to be made starting 
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Figure 3. - Typical infrared radiation oscillogram. 



at about 50 microseconds after 
passage of the shock.' 

It will he noted in figure 3 
that the gas behind the reflected 
shock in this particular run ab- 
sorbed almost all of the radi- 
ation from the glower source. In 
many runs of this sort, the por- 
tion of the absorption trace 
after the reflected shock was 
simultaneously displayed on an- 
other oscilloscope at much higher 
gain, so that the absorptivity 
could be accurately determined. 


Test Gas 

A commercially prepared 
argon- carbon dioxide mixture was 
used without further treatment. 

It analyzed 9. 9 percent carbon 
dioxide by volume; the rest was 
argon. A few tests were made 
with this gas diluted to 6.7 per- 
cent carbon dioxide in argon and 
with a commercially prepared 
argon- carbon dioxide mixture that 
analyzed at 1. 1 percent carbon 
dioxide and the balance, argon. 


RESULTS AND DISCUSSION 


Figure 4. - Comparison of measured and calcu- 
lated temperatures behind incident shocks in 
carbon dioxide-argon. 


The results of determina- 
tions of gas temperatures behind 
incident shocks are presented in 
figure 4, where they are compared 
with temperatures calculated by 

simple shock theory from measured shock velocities. The IMRA temperatures were 
actually measured at several 10-microsecond intervals for each experiment, 
starting at about 50 -microsecond laboratory time when the signals first became 
readable. Each point in figure 4 was obtained by extrapolating such data back 
to time zero, the instant when the shock passed the center of the windows. 

This was done to eliminate effects of shock velocity attenuation on the temper- 
ature. 

The agreement is good, and as expected, does not seem to be affected by 
the carbon dioxide concentration over the range studied. Despite this general 
agreement between measured and calculated incident temperatures, there is 
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nevertheless a good deal of scatter 
evident in figure 4. Part of this 
scatter is traceable to contamina- 
tion of the shock-tube windows, hut 
most of it is due to a drifting type 
of instability of the infrared de- 
tector. Prom run to run, this in- 
stability produced uncertainties as 

large as 10 percent in 1L . Equa- 

A, g 

tion (l) shows that this error re- 
sults in an uncertainty of only about 
30° K at temperatures near 1000° K, 
but the error grows to nearly 300° K 
at temperatures near 4000° K. 


Inasmuch as the main purpose of 
the work was to measure reflected 
shock temperatures at levels above 
2000° K, it was desirable to reduce 
the expected scatter in the results 
as much as possible. This was done 
by using the incident shock wave as an internal standard for each run. The 
temperature calculated from the measured shock velocity at the window position, 
and the measured emissivity of the gas behind the incident wave, were both as- 
sumed to be correct. The value of spectral radiance g required to satisfy 

equation (l) was then calculated. This value, and the measured voltages V 
(incident shock) and V Q , were inserted in equation (2) and used to calculate 
a lumped instrument factor, (l/R)(l/t v )N^ Q . 

In each run, then, this individually determined factor was used to convert 
measured voltages into a series of reflected shock temperatures, starting at 
about 50 microseconds behind the shock and determined at approximately 
10-microsecond intervals. These data were extrapolated to time zero, the in- 
stant the reflected shock passed the center of the windows. The resulting 
measured temperatures are plotted as ordinates in figure 5; the abscissas are 
the corresponding reflected temperatures graphically (ref. 10) calculated from 
the velocity of the incident shock as it passed the windows. 

Clearly, this means of obtaining calculated reflected shock temperatures 
for comparison with measured values is rather arbitrary. Its most obvious 
shortcoming is that it does not allow for the attenuation of the incident wave. 
However, it is important to see how well the temperatures calculated in this 
simple way agree with the measured values. Figure 5 shows that they agree 
quite well, although there are trends in the data. The points obtained with 
the end wall at its most remote position, 179 millimeters from the windows, 
tend to be high, and those obtained with the end wall 78 millimeters away tend 
to be low. This behavior is not inconsistent with that noted in reference 11, 
which reported that the reflected shock pressure behaved as if the shock decel- 
erated and then accelerated as it receded from the end wall. Despite these 
trends, the agreement between measured and simply- calculated temperatures is 
very gratifying. Of the 27 measurements shown in figure 5, 14 are within 



Figure 5. - Comparison of measured and calcu- 
lated temperatures behind reflected shocks in 
carbon dioxide-argon. 
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Figure 6. - Experimental and calculated temperatures behind 
shock waves in carbon dioxide-argon. = 51.0 torr ; in- 
cident Mach number, 3.64. 


100° K of the line, and only 4 miss it by more than 200° K. It is, however, 
believed that deviations larger than 100° K are real and that they reflect ac- 
tual departures of the shocks from ideal one-dimensional behavior. 


The two extreme points in figure 5 represent data from oscilloscope records 
in which the reflected portions could be accurately measured but the incident 
portions could not. These two points were therefore calculated with the origi- 
nal calibration factors rather than an individually determined factor such as 
was used for each other point. They are included to extend the temperature 
range of the reflected data and to indicate that the individual calibration cor- 
rections were usually small. 


An example of the detail which can be obtained by the IMRA method is given 
in figure 6, where the temperature history of a shock is followed at 10- to 20- 
microsecond intervals for 700 microseconds after its passage. Measurement of 
the incident shock wave velocity U at four successive places just before ar- 
rival at the test window showed the velocity to be attenuating at the rate of 
1.5x10“^ (millimeters per microsecond) per millimeter. Therefore, 
dU/U = 1.5x10"^ for a 1-microsecond interval. Strong shock theory predicts the 
relation 2dU/U = dT/T = dP/P. Therefore, in the first 200 microseconds after 
the passage of the incident shock, the temperature would be expected to rise 
from 1310° K (calculated from shock velocity) to about 1388° K and the pres- 
sure P from 1.76 to about 1.87 atmospheres. The IMRA measurements showed this 
rise approximately, as indicated by their conformity to the incident theory line 
in figure 6. The pressure rose to 1.90 atmospheres, or dP/P = 3.8x10"^ this 
also shows good agreement with the prediction. 


Other observers (ref. 2) have noted a rise in pressure behind reflected 
shocks. This pressure rise, converted into an isentropic temperature change 

by the equation T = T 0 (P/P Q )^ r w h ere y is the ratio of specific 
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Figure 7. - Experimental and calculated temperatures behind 
shock waves in carbon dioxide-argon. Pj = 51.0 t o r r ; in- 
cident Mach n umber, 4. 09. 


heats and T Q is the reflected shock temperature calculated from incident 
shock velocity, corresponds very veil vith the infrared pyrometer results of 
figure 6. 

A second example of detailed analysis is presented in figure 7. Here the 
temperature following the passage of the reflected shock tends to rise, because 
of isentropic compression, and tends to decrease because of dissociation accord- 
ing to the equation 


C0 2 -> CO + i 0 2 , AH = 65.8 kilocalories 

The reaction has atomic oxygen as an initial product. However, the recombina- 
tion rate of atomic oxygen is so high at these conditions that this step in the 
reaction mechanism can be ignored. The following were used to make theoretical 
temperature calculations: (l) 2.9x10^ cubic centimeters mole”^ second^ for 

the dissociation rate of carbon dioxide, calculated for 3220° K from the expres- 
sion given in reference 12, and for the subsequent temperatures after the re- 
flected shock-wave passage, (2) the isentropic temperature corrections from the 
measured pressures, and (3) negligible back reaction of the carbon monoxide and 
oxygen to form carbon dioxide. As with the simpler case of figure 6, the agree- 
ment of theory and results is gratifying. In this case, figure 7, the attenua- 
tion comparisons behind the incident shock are dU/U = 2.6x10“^, 

dP/P = 4.3X10" 4 , and dT/T = 3.4xl0" 4 . 
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CONCLUSIONS 


This study of gas temperatures measured in a shock tube hy infrared mono- 
chromatic radiation pyrometry has led to the following conclusions: 

1. Infrared pyrometry of a carbon dioxide-argon gas mixture yielded inci- 
dent shock temperatures in agreement with those calculated from the incident 
shock speed hy common shock theory, from 1100° to 3300° K. 

2. A comparison of reflected shock temperatures measured 27, 78, and 
179 millimeters after reflection with temperatures calculated from one- 
dimensional shock theory shows the latter are generally good for short dis- 
tances after reflection. 

3. Attenuation of the speed of an incident shock wave was accompanied hy 
changes in the gas pressure and temperature in the period following the passage 
of the shock. The relations between these changes were in general agreement 
with strong shock theory. 

4. The change in gas temperature with time after a reflected shock was 
adequately calculated from the dissociation rate of carbon dioxide and the pres- 
sure history of the gas mixture. 


Lewis Research Center, 

National Aeronautics and Space Administration, 
Cleveland, Ohio, June 4, 1965. 
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